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A simple scheme for conditional generation of nonclas- 
sical light with sub-Poissonian photon-number statistics is 
proposed. The method utilizes entanglement of signal and 
idler modes in two-mode squeezed vacuum state generated 
in optical parametric amplifier. A quadrature component of 
the idler mode is measured in balanced homodyne detector 
and only those experimental runs where the absolute value 
of the measured quadrature is higher than certain thresh- 
old are accepted. If the threshold is large enough then the 
conditional output state of signal mode exhibits reduction of 
photon-number fluctuations below the coherent-state level. 

PACS number(s): 42.50.Dv 



I. INTRODUCTION 

The possibility of a conditional generation of nonclas- 
sical states of light has been intensively investigated in 
recent years. The conditional quantum-state preparation 
schemes benefit from quantum entanglement between sig- 
nal mode and an ancilla system (e.g. two output ports 
of a beam-splitter H , signal and idler modes in paramet- 
ric down-conversion or cavity mode and an atom 
in cavity QED [Q). Suppose we measure some observ- 
able of the ancilla. The collapse of the ancilla state, 
caused by the measurement, influences, due to entan- 
glement, the state of the signal mode. In the conditional 
quantum-state preparation schemes we accept only those 
experimental runs where the required measurement out- 
come (or, more generally, a sequence of measurement out- 
comes) is observed, and we reject all unsuccessful trials. 

Based on this general strategy, schemes for conditional 
generation of Fock states Q , Schrodinger cat states |||7]] 
and arbitrary superpositions of Fock states in a cavity 
mode have been suggested. In particular, photon- 
number detection on idler mode of two-mode squeezed 
vacuum allows one to conditionally generate a highly 
nonclassical state of signal mode whose Wigner function 
may be negative near the origin of the phase space |P,^0[ . 
Alternatively, one can also use the detected intensity of 
idler beam as a negative feedback which modulates the 
intensity of pumping beam pj , or as a trigger of an op- 
tical shutter placed in path of pump or signal beams [ p"2| . 
All these techniques can provide a sub-Poissonian light 
in signal beam. 

A direct detection of the number of photons of the idler 
mode is not the only possibility here. A quadrature- 
measurement via balanced homodyne detection would 



prepare the signal in squeezed coherent state H||. Sim- 
ilarly, a heterodyne detection on idler leads to coherent 
state of the signal Unfortunately, the coherent am- 
plitudes are random, because they are proportional to 
the value of the detected quadrature of the idler mode. 
This problem may be circumvented by the feedforward 
technique where a suitable displacement of the signal is 
performed thereby compensating the undesired effects of 
the measurement on the idler mode Q . 

In this paper we propose an alternative simple and 
new way of conditional generation of nonclassical light 
in the signal mode via homodyne detection of the idler 
mode. In the suggested scheme, the experimental run is 
accepted only when the measured quadrature falls within 
a certain chosen window of quadratures. This guarantees 
that the conditional state is prepared with finite proba- 
bility, controlled by the size and position of the quadra- 
ture window. For the sake of concretness, we shall as- 
sume that the signal state is successfully generated when- 
ever the absolute value of the detected idler quadrature 
is higher than certain threshold xo. We shall see that 
the conditional state of signal can exhibit sub-Poissonian 
photon-number statistics. The amount of suppression 
of photon- number fluctuations can be controlled by xq. 
There is a trade-off between preparation-probability and 
the amount of photon-number fluctuations of the gener- 
ated state. Strong suppression of the fluctuations can 
be achieved at the expense of a small success rate and 
vice versa. We address the influence of imperfect detec- 
tion and show that the nonclassical light can be prepared 
only if the total efficiency of homodyning is higher than 
certain threshold. 

The paper is organized as follows. In Sec. II we provide 
a description of the suggested experimental setup. The 
quantum statistical properties of the conditional state 
of signal mode are studied in Sec. III. The influence of 
imperfect homodyne detection is discussed in Sec. IV. 
Finally, Sec. V contains conclusions. 

II. CONDITIONAL GENERATION 

Consider the experimental setup shown in Fig. 1. The 
two- mode squeezed vacuum state \ip) prepared in the 
nondegenerate parametric amplifier (NOPA) reads (l3| 
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where a s , <Zj (o|, at) are annihilation (creation) operators 
of signal and idler modes, respectively, \ri) s , \n)i are Fock 
states of the signal and idler modes, r is (real) squeez- 
ing constant of NOPA and A = tanh 2 r. We note that 
the currently available parametric amplifiers can produce 
strongly squeezed light. For example, Aytiir and Kumar 
Jl4| reported a parametric gain cosh 2 r f» 10, which cor- 
responds to A w 0.9. 

As illustrated in Fig. 1, the idler mode of two- mode 
squeezed vacuum is fed to the balanced homodyne detec- 
tor. The measurement on idler triggers the observation 
of the signal mode, which is performed only if the de- 
tected idler quadrature falls within the chosen window 
of quadratures X. In this way we conditionally generate 
the output state of the signal mode. 

Since it is experimentally challenging to maintain a 
fixed phase of the local oscillator, we consider the ho- 
modyne detection with randomized phases. Positive 
operator-valued measure (POVM) describing this type 
of quantum measurements reads 
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where \x; 4>)i is an eigenstate of the quadrature operator 

Xi,4, = {Sue-** + a\e^)/V2, 



(3) 



The POVM corresponding to the quadrature window X 
can be expressed as 
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The (normalized) reduced density matrix p s of the signal 
can be obtained by tracing over idler mode 
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where I s is unit operator on Hilbert space of signal mode. 

Due to the phase randomized homodyning, density ma- 
trix of the conditional state of the signal mode is diagonal 
in Fock basis, 
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On inserting the POVM (@) into Eq. (|), we get formula 
for the photon- number distribution p n , 
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where 



C(A, X) = (ip\l s <g> ftilV) = (1 - A) X n q n (8) 
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FIG. 1. Scheme of conditional generation of sub-Poisso- 
nian light in signal mode from two-mode squeezed vacuum 
prepared in nondegenerate parametric amplifier (NOPA). A 
quadrature of idler mode is measured in balanced homodyne 
detector, which consists of a beam splitter BS, strong coherent 
local oscillator LO, two photodetectors PD and a subtracter 
of the two signals. Detected quadrature falling within the 
chosen window X triggers the observation of the signal mode. 



is a normalization coefficient and 
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Here H n (x) denotes Hermite polynomial of variable x. 



III. PHOTON-NUMBER STATISTICS OF THE 
CONDITIONAL STATE 

So far we have assumed an arbitrary X. In what fol- 
lows, we analyze the case where the conditional state of 
signal mode is successfully generated if the absolute value 
of the idler quadrature is higher than certain lower bound 
xq. In this case we can write, 



Hl{x)e~ x dx. 



(10) 



The normalization factor C(A,a;o) represents the proba- 
bility of conditional generation of the output state (|5|), i.e. 
the probability that the absolute value of the measured 
quadrature of the idler would be larger than xq. The 
idler mode is in thermal (chaotic) state and its quadra- 
ture component exhibits Gaussian distribution with 
zero mean and variance (1+A)/[2(1 — A)]. The probability 
that | a; | > xq can be thus expressed as 



exp —x 



After a straightforward integration we arrive at 
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where the error function is defined as 
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FIG. 2. Mean photon number (h) and Mandel Q-factor 
versus xo for three different squeezing parameters A = 0.05 
(solid line), A = 0.1 (dashed line) and A = 0.2 (dot-dashed 
line) . 



The function C(A, xo) contains a complete information 
about the quantum-statistical properties of the state (|^). 
From the definition (Q) we can deduce that C(A, Xq) is 
a generating function of moments of the photon number 
distribution (Q), 
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Also the normally ordered moments can be determined 
from C(A, Xq), 
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In particular, we obtain expressions for the mean photon 
number 
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and second factorial moment 
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The conditional output states p s can be highly nonclas- 
sical, exhibiting sub-Poissonian statistics. The suppres- 
sion of photon-number fluctuations can be conveniently 
characterized by Mandel Q-factor p5[ 
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FIG. 3. (a) The probability of conditional generation of 
state with Q — (solid line), Q = —0.05 (dashed line) 
Q — —0.1 (dot-dashed line), and Q = —0.2 (short-dashed 
line). Figure (b) shows the corresponding values of xo- 



It holds that Q > — 1 and the equality is achieved for 
Fock state. The light is sub-Poissonian when the photon- 
number variance ((An) 2 ) is less than the mean number 
of photons (h). This is indicated by negative value of 
Q. The statistics are Poissonian when Q = and super- 
Poissonian if Q > 0. 

The dependence of the mean photon number ( |l6| ) and 
the Q-factor (||) on the lower bound £o is plotted in 
Fig. 2 for three different values of A. We can see that 
(n) monotonically grows with xq while Q monotonically 
decreases and eventually becomes negative. The fig- 
ure clearly indicates that the states p s can possess sub- 
Poissonian statistics. However, we should be concerned 
with the probability C(X,xq) of generation of such non- 
classical states. 

Figure 3a shows dependence of this probability on A 
for four different Q-factors. The corresponding values of 
xq required for reaching given Q are displayed in Fig. 
3b. Consider first the boundary case when the photon- 
number fluctuations are suppressed to the level of Pois- 
sonian distribution, hence Q = 0. The probability of 
generation can be almost 55% for weakly squeezed states 
(A <C 1) and decreases with growing A. In agreement with 
Fig. 2b, the xq required for achieving Q — increases 
with A. It is interesting that there is a lower bound on 
xo as A — * 0. In order to suppress the fluctuations to 
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FIG. 4. Photon-number distribution of the conditional 
state p s for A = 0.25 and (a) xq = 0, (b) a;o = 1, (c) xq = 2, 
and (d) xo = 3. 



and rewrite Eq. (|h]) as 
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After some algebra we find that 
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On performing the differentiation we finally arrive at 
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This expression looks particularly simply when rewritten 
as a recurrence relation, 



Poissonian level, we need finite Xq even for very weak 
squeezing. This minimum .To im i n can be determined if 
we consider the formulas (pf), (17) and ( |l8| ) in the limit 
A — > 0. It follows that Xo,min is a root of a transcendent 
equation 
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whose numerical solution provides XQ^nn — 0.4248. 

Assume now that we want to generate a nonclassical 
state with negative Q. Figure 3 reveals that for each 
Q < there is certain optimal value of A which allows us 
to generate this state with highest probability. The peak 
in probability roughly corresponds to the minimum of xq 
but those two extrema do not coincide. 

Having investigated the photon-number variance we 
turn our attention to the photon-number distribution it- 
self. The coefficients q n can be calculated from C(X,xq) 
as follows, 
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This formula, however, is not very convenient for practi- 
cal calculations. Therefore we present a different deriva- 
tion of q n which leads to formulas suitable for numer- 
ical processing. In order to perform the integration in 
Eq. ([To]), we invoke the generating function for Hermite 
polynomials, 
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starting with qo — 1 — erf(xo). 

Figure 4 illustrates how the shape of the photon num- 
ber distribution of p s changes with xq- When xq = 
and all quadratures are accepted, the signal mode is 
in thermal state with mean number of chaotic photons 
(h) = A/(l — A) and positive Q-factor (Q a = 0.333), 
see Fig. 4a. This initially super-Poissonian distribu- 
tion becomes narrower when the lower bound of accepted 
quadratures xq is increased and the statistics shown in 
Figs 4b-d are sub-Poissonian, as indicated by negative Q- 
factors (Q b = -0.026, Q c = -0.216, and Q d = -0.297). 
The selection conditioned by quadrature measurements 
suppresses the contribution of vacuum and low Fock 
states and p n becomes peaked about n which grows with 

Xq. 

These changes in the form of the photon-number dis- 
tribution are reflected in the shape of quasidistributions, 
which are phase space representations of the quantum 
state. The quasidistributions for p s are axially symmet- 
ric around the origin of phase space because the density 
matrix p s is diagonal in Fock basis and does not carry 
any information about phase. The Husimi quasidistribu- 
tion W^i(a) = ir^ 1 (a\p s \a) corresponding to states from 
Figs. 4a and 4c is plotted in Fig. 5. The quasidistri- 
bution of thermal (chaotic) state shown in Fig. 5a has 
Gaussian shape and is peaked at the origin of the phase 
space. On the other hand, depicted in Fig. 5b has a 
clear dip in the origin and is peaked at a nonzero radius. 

Although the conditional state p s can be nonclassical, 
its Wigner function W s (x s ,p s ) is always positive. This 
should be contrasted with conditional generation based 
on photodetection on the idler, which may lead to states 
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FIG. 5. The Husimi quasidistributions of the state p s for 
A = 0.25 and (a) x = 0, (b) x = 2. 



with negative Wigner function |9[|1C[]. The positivity of 
W s can be proved by recalling that the trace of two op- 
erators can be evaluated by integrating the product of 
their Wigner functions over the whole phase space. The 
Wigner representations of the two-mode squeezed vac- 
uum (|l|) and POVM (H) are both positive. Taking into 
account the definition (|]) of p s , we immediately conclude 
that the Wigner function W s of the state p s must be 
positive. In this context we can say that direct pho- 
todetection by avalanche photodiode is, in certain sense, 
more non-classical than balanced homodyning because 
the Wigner function of the operator II = I — |0)(0| is 
negative in some regions of phase space. 



IV. IMPERFECT HOMODYNE DETECTION 

So far, we have assumed an ideal homodyne detector 
with unit efficiency. Let us now analyze the influence of 
imperfect detection with efficiency 77 < 1 on the proper- 
ties of conditional state of signal mode. We can expect 
that the nonclassical properties of p s will diminish with 
decreasing 77 but we shall see that the suppression of non- 
classicality is not severe. 

A realistic homodyne detector can be modelled by an 
ideal homodyne detector whose signal input is preceded 
by a beam splitter, where the signal mode is mixed with 



an auxiliary mode which is in chaotic state. The detected 
quadrature can be thus expressed as |16| 



(25) 



In our case the quadrature Xi tVt j, has again a Gaussian 
distribution with zero mean and variance which does not 
depend on <^>, 
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where n is mean number of chaotic photons in the auxil- 
iary mode. Making use of the formula (|2^) we can derive 
the probability that the absolute value of the measured 
idler quadrature is higher than the threshold x , 



C (A, xq , 77, n) — 1 — erf 
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In what follows we shall assume that the auxiliary 
mode is in vacuum state, hence n = and Eq. (|2T 
simplifies to 



C(X,x ,rj) = 1 - erf a; 
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In the limit of ideal detector, 77 = 1, this expression re- 
duces to Eq. (H). We may utilize the formulas jl4|), 
15|) and calculate the moments of the photon number 
stribution (Q) . After some algebra we obtain 
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The POVM describing homodyne detection with real- 
istic detector can be expressed as 



n, 



where Ilj(x^) reads 
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The coefficients q n in the photon-number distribution (j?]) 
may be again calculated in a similar way as in the case of 
ideal detector. After a long but straightforward algebra 
we arrive at recurrence relation for q n , 
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where we start from qg = 1 — erf(xo). 

The detrimental effect of rj < 1 on the suppression of 
photon-number fluctuations can be deduced from Fig. 6. 
In Fig. 6a we plot the probability C(X,xq,t]) of gener- 
ation of the state with Poissonian fluctuations (Q = 0). 
We can see that the probability reduces with decreas- 
ing rj. Our numerical calculations reveal that we can- 
not achieve Q = for any A when rj < 0.5. Figure 6b 
shows the probability of preparation of a state with sub- 
Poissonian statistics (Q — —0.05). As rj diminishes, the 
maximum probability of generation decreases and the op- 
timal A increases. 

We can provide a simple argument explaining why non- 
classical light cannot be generated when 77 < 1/2. To 
this point let us assume that the homodyne detector in 
Fig. 1 is replaced with heterodyne detector measuring 
projections to coherent states, IT (a) = ir~~ 1 \a)i(a\. Sup- 
pose that the output signal is accepted only when the 
absolute value of the real part of measured a is larger 
than ao > 0. Since the idler is in chaotic state, its Q- 
function is Gaussian and we may express the probability 
that |Re(a)| > ao as 



C(A,a ) 



2(1 -A) 



C(A, a ) = 1 - erf (aoVl^A). 



(34) 



We can see that C*(A,ao) = C(A,xo = ao,?7 = 0.5) 
hence the photon-number statistics of conditional sig- 
nal state generated via heterodyne detection on idler 
coincides with statistics of conditional signal state pre- 
pared by means of homodyne detection on idler with effi- 
ciency 0.5. After the detection of IT (a) the signal mode 
is projected to coherent state with complex amplitude 
a s = VAa* |^| . Thus the conditional state prepared by 
heterodyne detection on idler is a stochastic mixture of 
coherent states and cannot exhibit any nonclassical fea- 
tures. Assume now that the auxiliary mode in Eq. (|2T 
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FIG. 6. The probability of generation of state with Man- 
del Q-factor (a) Q = and (b) Q = —0.05 is plotted 
for four different detection efficiencies r\ = 0.9 (solid line), 
7j — 0.8 (dashed line), rj — 0.7 (dot-dashed line) and rj = 0.6 
(short-dashed line). 




is in thermal state. When we look at the formula (|2J 
we find that C(A, Xq, rj, n) has the same form as C(A, ao) 
in Eq. (||) when 



V = Vth 



1 + 2n 

2 + 2n' 
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If rj < ?/th, then the balanced homodyning is formally 
equivalent to heterodyning with detection efficiency lower 
than unity. Thus we conclude that the sub-Poissionian 
light can be conditionally generated only when rj > rj t h 
holds. The threshold efficiency 77th is lowest when the 
auxiliary mode is in vacuum state and rj t h = 0-5- The 
threshold 77 t h increases with n and approaches 1 for large 
n. 

In the experiments, the efficiency is reduced mainly 
due to only partial overlap between the measured mode 
and local-oscillator mode. If the overlap is good, then ?y 
can be very high. For example, in the single-mode opti- 
cal homodyne tomography ]17j , the efficiency 82% ± 5% 
was achieved, which would be fully sufficient for our pur- 
poses. However, we deal with two-mode system and we 
need to extract the signal mode which corresponds to de- 
tected idler. In the recent experimental realization of the 
two-mode homodyne tomography [ p"8[ , the estimated effi- 
ciency was only 35%, which is below our threshold. Nev- 
ertheless, even such data may be used to demonstrate the 
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suppression of noise via conditional generation, although 
the photon number fluctuations could not be reduced be- 
low the Poissonian level. 



V. CONCLUSIONS 

We have proposed a simple scheme for conditional 
generation of sub-Poissonian light from the two-mode 
squeezed vacuum. The scheme involves homodyne detec- 
tion on the idler mode and the experimental run is ac- 
cepted only if the absolute value of the detected quadra- 
ture is higher than certain threshold. As the threshold is 
increased, the photon-number fluctuations of condition- 
ally generated state of signal mode are gradually sup- 
pressed but the probability of generation decreases. If 
we wish to prepare a nonclassical state with a given neg- 
ative Mandel Q-f actor, then we can choose an optimal 
squeezing parameter A leading to maximum generation 
probability. Nonclassical light can be prepared only when 
the efficiency of homodyne detection rj is higher than 
Vth > 1/2. Although the threshold 1/2 was not achieved 
in the recent experiment on two-mode homodyne tomog- 
raphy 18 1, the condition jj > 77th is by no means severe 
and the present method may become experimentally fea- 
sible in a near future. 
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